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ABSTRACT 

The  thesis  provides  two  subroutines  for  the  Interim 
Battle  Group  Tactical  Trainer  (IBGIT).  The  first  subroutine 
is  a  single  sensor  tracking  model  using  the  Kalman  filter. 
This  subroutine  is  part  of  the  Passive  Sonar  Model.  The 
second  subroutine  is  a  multiple  sensor  tracking  model  using 
the  Kalman  filter  to  correlate  all  sonar  contacts  on  a 
specific  target.  This  subroutine  is  a  separate  entity  and 
can  be  turned  on  or  off  at  simulation  initiation  as  required 
by  training  objectives. 
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THESIS  DISCLAIMER 


The  reader  is  cautioned  that  computer  programs  developed 
in  this  research  may  not  have  been  exercised  for  all  cases 
of  interest.  While  every  effort  has  been  made,  within  the 
time  available,  to  ensure  that  the  programs  are  free  of 
computational  and  logic  errors,  they  cannot  be  considered 
validated.  Any  application  of  these  programs  without  addi¬ 
tional  verification  is  at  the  risk  of  the  user. 
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I.  INTRODUCTION 


The  Naval  Ocean  Systems  Center  has  developed  the  Interim 
Battle  Group  Tactical  Trainer/Computer  Support  Facility 
(IBGTT/CSF)  as  a  computer-based  tactical  simulation  system 
to  provide  a  training  device  for  senior  naval  officers  to 
practice  tactical  decision  making  until  such  time  as  the 
Enhanced  Naval  Warfare  Gaming  System  becomes  available.  The 
trainer  is  intended  to  provide  inte~active,  multithreat, 
multiplatform  operational  situations  in  a  simulated  yet 
realistic  operational  environment  so  that  selected  officers 
can  study,  practice,  and  be  evaluated  in  force-level 
tactical  decision  making. 

The  IBGTT  training  capability  is  implemented  as  a  real¬ 
time,  man-interactive,  computer-aided  (discrete  event,  time 
step)  simulation  of  the  naval  warfare  environment.  In  opera¬ 
tion,  the  IBGTT  supports  a  two-sided  (BLUE  vs.  ORAm'GEji 
interactive  scenario  in  which  opposing  sides  can  define, 
structure,  and  dynamically  control  forces  ranging  in  size 
from  one  or  more  battle  groups  an  .  associated  aircraft,  down 
to  a  single  air  or  surface  unit.  Force  elements  and  their 
associated  characteristics,  sensors,  weapons,  and  communica¬ 
tion  systems  may  be  derived  iron  real,  proposed,  or  concep¬ 
tualized  units  or  systems. 

The  utilization  of  IBGTT  involves  the  use  of  the  four 
major  Naval  Warfare  Interactive  Simulation  System  (NWISS) 
processes;  BUILD,  FORCE,  WARGAME,  and  POST-GAME  ANALYSIS. 
The  BOILD  process  is  a  stand  alone  interactive  program  used 
to  create  and  maintain  platform,  sensor,  communication,  and 
weapon  characteristics  in  the  IBGTT  Characteristic  Data 
Base.  The  FORCE  process  is  a  stand  alone  interactive  program 
used  to  create  and  maintain  an  exercise  scenario  using  the 


data  base  created  by  the  BUILD  process.  The  WARGAME  process 
is  an  interactive  program  used  to  accept  and  execute  user 
orders;  control  platform  motion,  detections,  and  communica¬ 
tions;  determine  engagement  and  other  event  outcomes;  and 
display  status  information  and  tactical  situations.  The 
POST-GAME  ANALYSIS  process  analyzes  and  lists  critical  data 
recorded  during  the  exercise;  supports  exercise  reconstruc¬ 
tion;  and  tentatively  evaluates  some  Measures  of 
Effectiveness. 

A  global  data  area  in  N1TISS,  called  the  blackboard,  is 
shared  by  all  the  major  modules  functioning  during  the  exer¬ 
cise.  It  is  the  area  where  these  modules  interface  with  each 
other  through  the  application  of  uniform  naming  conventions 
and  the  efficient  use  of  memory.  The  blackboard  is  essen¬ 
tially  comprised  of  numerous  tables  and  subtables.  Each 
table  is  assigned  specific  pointers  while  each  subtable  is 
assigned  specific  indices.  The  tables  and  subtables  contain 
the  fields  (data)  that  are  unique  to  that  particular  table 
or  suitable.  Each  data  item  in  the  blackboard  is  referred  to 
as  a  field  and  includes  both  whole  words  and  specific  bits. 
The  field  names  are  structured  to  provide  tne  identity  of 
the  associated  pointers  and  indices  as  well  as  the  data  type 
in  the  field. 

Effective  training  at  this  level  requires  models  of 
naval  warfare  interactions  which  provide  realistic  results 
based  on  an  emulation  of  the  actual  warfare  system.  NWISS 
uses  a  wide  variety  of  models  to  simulate  the  behavior  of 
platforms,  weapons,  sensors,  and  communication  systems. 
However,  many  of  these  models  do  not  emulate  the  actual 
warfare  system  nor  do  they  provide  realistic  results. 
Therefore,  it  is  necessary  to  improve  or  replace  these  defi¬ 
cient  models  in  order  to  obtain  effective  training  and  meet 
the  objectives  for  which  IBGTT  was  designed. 
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This  thesis  will  address  two  models  in  particular.  The 
first  model  is  the  Target  Motion  Analysis  (TMA)  Model  which 
processes  passive  sonar  contacts.  The  current  model  will  be 
examined,  followed  by  a  presentation  of  a  Kalman  filter 
improvement  to  the  model.  Secondly,  the  Sonar  Correlation 
Model  will  be  examined,  followed  by  a  presentation  of  a 
Kalman  filter  to  replace  the  current  model. 

The  reader  should  be  advised  that  the  TMA  Kalman  Filter 
Model  is  in  actuality  equivalent  to  the  Sonar  Correlation 
Kalman  Filter  Model,  as  will  become  evident  in  the  presenta¬ 
tion  of  the  two  Kalman  Filter  Models.  This  thesis  further 
presupposes  that  the  reader  is  familiar  with  the  Kalman 
filter. 


II-  SINGLE  SENSOR  TRACKING 


A.  TARGET  MOTION  ANALYSIS  (TMA)  MODEL 

The  current  model  will  monitor  the  number  of  game 
minutes  for  which  passive  contact  has  been  held  on  each 
target  by  each  detecting  passive  sonar  (i.e.,  submarine, 
surface  sonar,  towed  array,  or  sonobuoy)  .  When  this  time 
exceeds  the  TMA  time  defined  by  the  user  at  simulation 
initiation  a  target  motion  analysis  report  will  be  displayed 
or  the  Passive  Sonar  Status  Board.  The  TMA  range,  course, 
and  speed  displayed  are  derived  as  follows: 

range"]  range" 

TMA  course  =  Actual  target  course  ±  FACTOR 
speed  speed 

The  FACTOR  is  the  product  of  a  random  number  drawn  from  a 
uniform  distribution  and  a  derived  parameter. 

These  parameters,  which  are  indicated  in  Table  I,  cause 
increasingly  accurate  solutions  to  be  developed  as  Signal 
Excess  (SE)  increases  and  as  the  target’s  true  bearing 
changes  (*B=delta  B)  from  the  true  bearing  of  its  initial 
detection.  This  latter  factor  simulates  improved  solutions 
derived  from  highec  bearing  rate  targets  and  longer  tracking 
times.  The  solution  quality  displayed  will  be  selected  from 
Table  II  as  a  function  of  SE  and  *3. 

Once  a  TMA  solution  has  been  displayed,  only  the  range 
is  updated  on  the  display.  The  range  only  update  continues 
until  the  signal  excess  and/or  change  in  true  bearing  cause 
a  new  parameter  to  be  developed  from  the  table  (e.g.,  SE 
changes  from  -6  to  -5  or  »B  changes  from  5  to  6)  .  When  a  new 
parameter  is  selected,  a  new  TMA  range,  course,  and  quality 
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IT.  TEST  EESOLTS 


The  single  sensor  and  multiple  sensor  models  have  been 
tested  on  a  limited  basis.  Each  subroutine  has  been  tested 
independently  to  insure  that  the  models  perform  as  designed. 
However,  the  subroutines  have  rot  been  tested  for  integra¬ 
tion  into  the  overall  NWISS  program  or  other  subroutines. 

Due  to  time  constraints  and  computer  availability,  inte¬ 
gration  tests  were  not  possible.  The  added  blackboard  space 
reguirea  has  not  been  made  but  should  not  pose  any  problems. 
Inherent  with  any  new  subroutine  is  the  unforeseeable  affect 
it  may  have  on  unrelated  subroutines.  This  aspect  of  testing 
has  not  been  performed. 


The  track  quality  will  be  the  same  as  that  presented  in 
Chapter  Two. 


D.  C BADGES  TO  OTHEB  SUBROUTINES 

The  rational  FORTRAN  and  source  code  listing  for  this 
model  are  contained  in  Appendix  C  and  D,  respectively.  The 
changes  presented  in  Chapter  Two  for  Subroutine  WARCYC  are 
also  applicable  to  this  model.  The  only  other  change  that 
will  be  necessary  is  that  all  active  information  needs  to  be 
added  to  the  Remote  Table. 
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no  fix  or  track  is  associated  with  the  correlation.  As  a 
result,  each  game  minute  two  bearing  lines  are  displayed 
(they  may  not  be  the  same  two  from  the  previous  minute  nor 
necessarily  an  improvement)  that  jump  around  the  screen, 
providing  no  useful  information  to  the  user. 

The  Kalman  Filter  Model  permits  all  the  information 
available,  both  active  and  passive,  to  be  correlated  on  a 
specific  target  and  be  displayed  as  a  fix  with  an  updated 
track.  In  addition,  the  track  quality  associated  with  the 
fix  provides  the  user  with  the  added  information  about  the 
relative  size  of  the  AO?.  The  on  and  off  ability  of  the 
Sonar  Correlation  Model  will  be  incorporated  into  the  Kalman 
Filter  Model  for  the  previously  stated  reasons  concerning 
its  flexibility. 

C.  KALMAN  FILTEB  MODEL 

As  mentioned  in  the  introduction,  the  single  sensor 
tracking  model  is  in  actuality  equivalent  to  the  multiple 
sensor  tracking  model.  The  difference  being  the  scope  of  the 
information  being  processed.  The  single  sensor  model  is  a 
subset  of  the  Passive  Sonar  Model;  whereas,  the  Multiple 
sensor  model  is  a  separate  entity  correlating  all  available 
information.  The  basics  of  the  models,  including  assumptions 
and  initial  conditions,  are  the  same.  Therefore,  only  the 
differences  from  the  single  sensor  tracking  model,  presented 
in  Chapter  Two,  will  be  discussed  here. 

The  model  will  handle  bearing  and  range  measurements  as 
well  as  bearing  only  measurements.  The  bearings  are  ±.5 
degrees  and  the  ranges  are  ±.5  nautical  miles.  For  the 
bearing  only  measurement,  the  K  matrix  will  be  the  same  as 
for  the  single  sensor  tracking  model.  For  the  bearing  and 
range  measurement,  the  following  H  matrix  will  be  utilized: 

H  =  f  -sin  (theta) /range  cos  (theta) /range  0  0 
I  cos  (theta)  sin  (theta)  0  0 


III.  HOITIPLB  SENSOR  TRACKING 


A.  SONAR  CORRELATION  MODEL 

The  Sonar  Correlation  Model  would  be  more  appropriately 
described  as  a  procedure  instead  of  a  model.  This  routine 
determines  the  correlation  of  bearings  between  two  detecting 
units  at  a  specific  target  and  performs  two  functions. 
First,  the  two  detecting  units  (with  an  intersection  angle 
of  at  least  60  degrees,  or  else  the  largest  available  angle) 
to  a  common  target  will  display  bearing  lines.  Second, 
multiple  targets  detected  within  a  certain  maximum  arc  will 
display  only  one  line;  and  will  set  the  composition  field 
set  (i.e.,  one,  few,  or  many).  All  other  passive  sonar  lines 
will  not  be  displayed. 

The  Sonar  Correlation  Model  can  be  turned  on  or  off. 
This  allows  for  the  flexibility  of  being  utilized  for  Battle 
Group  Commanders  and  their  staff  when  the  "big  picture"  is 
the  main  concern  and  not  the  individual  unit  prosecution; 
and  yet,  be  turned  off  when  the  trainer  is  being  utilized 
for  a  single  unit  or  group  of  units  practicing  coordinated 
operations. 

B.  KALMAN  FILTER  REPLACEMENT 

The  current  procedure  has  many  drawbacks.  Only  passive 
sonar  lines  are  considered;  active  information  is  displayed 
separately  and  is  not  correlated  with  the  passive 
informa tion. 

The  routine  searches  through  the  Remote  Table  until  it 
finds  two  bearing  lines  that  meet  the  60  degree  criteria. 
These  are  displayed  and  the  routine  ceases  to  search;  thus, 
not  necessarily  choosing  the  optimum  solution.  In  addition, 
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2.  Subroutine  LCLPSN 
- - *"* —  — — —  — — — — — 

Change  line  9  to  read: 

EQOI  VALENCE  {IBB,  FBB,  CBE,  IBBW,  I3BB, PBB) 

This  includes  the  blackboard  of  the  P  matrix. 
After  line  113  add: 

I_LCL$LASTTMATIME=  (IAND(ISHPI  (IBB (KPOINT_LCL 
*+1J  ,-0)  ,65535) 

This  time  is  used  in  Subroutine  LCLTMA. 

Change  line  262  to  read: 

IF  (.  NOT.  (I_ICL$OMN  IFL  AG.  NE.  1)  )  GOTO  23269 

This  removes  the  TMA  exceed  time  criteria. 

Change  line  264  to  read: 

LCLTMA (KPOINT_LCL , RO  BLAT, HOHLON ,LBEAR , 
*I_LCL$PMATBIX,I  LCL$ IASTTMATIME) 


•  If  the  2-sigma  axis  is  less  than  or  equal  to  500  yards, 
then  the  track  receives  a  guality  of  GOOD.  The  criteria 
of  500  yards  is  used  because  the  Engagement  Model 
employs  a  500  yard  kill  radius  for  a  torpedo. 

•  If  the  2-sigma  axis  is  greater  than  500  yards  but  less 
than  or  equal  to  1000  yards,  then  the  track  receives  a 
quality  of  FAIR.  The  criteria  of  1000  yards  was  chosen 
simply  because  it  is  twice  the  GOOD  criteria. 

•  If  the  2-sigma  axis  is  greater  than  1000  yards,  then 
the  track  receives  a  quality  of  POOR. 

E.  CHANGES  TO  OTHER  SUBROUTINES 

Implementation  of  this  subroutine  (LCLTMA,  rational 
FORTRAN  and  source  code  listing  are  contained  in  Appendix  A 
and  B,  respectively)  will  require  some  changes  to  other 
subroutines  and  additions  to  the  blackboard.  These  changes 
and/or  additions  include: 

1 .  Subroutine  WAECYC 

Each  target  needs  a  P(0)  matrix  in  the  blackboard; 
therefore,  add: 

REAL  PBB  (4,4) 

REAL  I_LCL$P MATRIX (4,4) , IJ3MT SPMATRIX  (4,4) 

DO  50000  J=1 , 4 

DO  50001  K=1 , 4 

IF( -NOT. (J.EQ. K)  )  GOTO  50001 

PBB  (J,K)  =  1000. 

GOTO  50000 
50001  PBB  (J,K)  =  0. 

50000  CONTINUE 

I_LCL$PMATEIX  =  PBB 
I  RMTJPMATRIX  =  PBB 
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Q0  is  the  initial  observation 


P<0)  = 


0  1000 


1000 


|_  0  0  0  1000  J 

Note:  1000  was  chosen  since  it  is  approximately  equal  to  ± 

one  convergence  zone  (32  nm)  and  z  32  knots. 

The  first  assumption  is  that 

E[W  (k)  *V  ( j)  Itranspose  ]  =  0  for  all  j  and  k. 

This  means  that  the  plant  noise  and  the  measurement  noise 
are  uncorrelated.  Secondly,  recall  that  the  assumption  is 
that  during  an  encounter,  the  target's  course  and  speed 
remain  constant. 

D.  TRACK  QUALITY 

In  order  to  reduce  the  number  of  changes  required  to  the 
overall  program,  the  TMA  quality  currently  used  from  Table 
II  will  be  utilized  tut  based  on  different  criteria.  This 
will  eliminate  the  need  to  change  the  blackboard;  and  more 
importantly,  will  not  change  the  Status  Tableau  seen  by  the 
player,  which  is  already  full. 

The  track  quality  is  based  on  the  semi-major  axis  of  the 
AOP.  This  is  determined  from  the  error  covariance  matrix  as 
follows: 

(semi-major  axis)  2  =  (pi  1*p22) /2  +  SQHT[  (pi  1-p22)  2/4+p122  ] 

where  pij  are  the  elements  of  the  P(t)  matrix.  A  2-sigma 
semi-major  axis  is  used  to  insure  a  probability  of  3.  3647. 
The  track  quality  is  then  deter  lined  as  follows: 


V  (t)  is  the  measirement  noise.  It  is  approximately 
N  (0 # E  (t) }  .  For  this  model,  all  bearings  are  ±  .5  degrees. 

K  (t)  is  the  Kalman  gain.  Ihe  update  from  X3hat  (t )  just 
before  the  measurement  to  X$hat  (t)  just  after  the  measure¬ 
ment  is  proportional  to  the  r^hcch;  the  Kalman  gaii*  is  the 
proportionality  constant. 


It  was  stated  earlier  that  the  measurements  are  assumed 
to  be  linearly  related  to  the  system  state.  Since  the 
measurement  is  in  polar  coordinates,  h(x)  is  in  fact 
nonlinear.  Therefore,  a  transformation  must  be  made  on  h  (x) 
to  give  a  linearly  related  H. 

In  this  model,  the  observation  will  be  made  from  a  plat¬ 
form  at  (u,w)  to  a  target  at  (x,y)  ,  where  x  is  north  and  y 
is  east.  So, 

h  (X)  =  theta  =  tan-*[  (y-w)  /  (x-u)  ] 
or, 

H  =  [  3h(X)/3x  3h(X)/3y  3h(X)/3x’  3h  (X)  /3y'  ] 

a  represents  the  partial  derivative 

evaluated  at  X  =  X$hat.  Hence, 

H  =  [  -sin  (theta) /range  cos  (theta)  /range  0  0] 

This  model  was  built  upon  two  initial  conditions  and  two 
important  assumptions.  The  initial  conditions  are: 

E[X  (0)  ]  =  X$hat  (0) 

and 


E[  (X  (0)  -XShat  (0 ) } *  [X  (0)  -X$hat  (0) }  ^transpose  ]=  P  (0) 

where 


X$hat(0)  =  [  32cos@0  32sin0e  0  0  ]$transpose 


P  is  the  error  covar;  ce  natrix. 

P(t)  =  E[  £X  (t)-X$h  (t)  }*  tX  (t) -X3hat  (t) }  $transpose  ] 

W  (t)  is  the  plant  noise.  It  describes  the  randomness  of 
the  system  as  it  moves  from  state  X(t)  to  X(t+1).  W  (t)  is 
approximately  N(0,Q{t)).  For  this  model,  Q  is  taken  to  be 
zero. 

Next,  a  new  fix  is  computed  based  on  an  observation. 
This  is  determined  from  the  measurement  model: 

Z  (t)  =  H(t)*X  (t)  ♦  V  (t) 

Thus,  mea  irement  is: 

Kalman  Gain:  K(t)  = 

P  (t)  *H$  transpose  (t)*[H(t)*P(t)  *H$ transpose  (t) 

+  H  { t)  ]-», 

State  Update: 

XShat  (t)  :=  X$hat  (t)  ♦  K  (t)  *[  z ( t)  -H  (t)  *X$hat  (t)  ], 

and 

Error  Covariance  Update: 

P(t)  :=  P  (t)  -K  (t)  *[  P  (t)  *H$transpose  (t)  ]$  transpose 

where 

:=  indicates  that  the  right  hand  side  is  computed  and 
replaces  the  value  on  the  left  hand  side  of  the  symbol. 

Z(t)  is  the  actual  measurement.  The  measurements  are 
assumed  to  be  linearly  related  to  the  system  state  X  ( t )  by 
the  observation  matrix  H  (t)  .  Note:  H  (t)  *X$hat  (t)  is  the 
predicted  outcome  of  the  measurement.  The  difference, 
Z (t) -H  (t) *X$hat  (t)  ,  is  the  measurement  residual  or  shock. 


problem  since  the  computed  AOE  updates  smoothly  with  the 
sensor  information. 


C.  KALHAH  FILTER  HODEL 

First  of  all,  it  is  assumed  that  during  an  encounter, 
the  target's  course  and  speed  remain  constant.  The  model 
updates  the  position  of  the  fix  since  the  last  observation 
base!  on  the  previous  estimate.  This  is  based  on  the  system 
model: 

X  (t)  =  PHI(t-1)*X  (t-1)+W(t-1) 

Thus,  movement  is: 

State  Extrapolation:  X$hat  (t)  =  PHI  (t-1)  *X$hat  (t- 1) 


Error  Covariance  Extrapolation: 

P(t)  =  PHI  (t-1)  *P  (t-1)*PHI$transpose  (t-1)  +  q  (t-1 ) 

where 

X$hat  is  the  estimated  state  vector.  It  is  assumed  to 
be  multivariate  normal  with  mean  zero. 

X (t)  =  [  x  (t)  y(t)  x'  y'  ]$transpose 
x  '=  x  velocity 
y'=  y  velocity 

PHI  is  the  transition  matrix.  It  describes  how  the 
state  vector  changes  from  X(t)  to  X  (t+1)  . 


PHI  =  1  0  «t  0 

0  1  0  •  t 

0  0  10 
0  0  0  1 


•t=delta  t 


B.  KALMAR  FILTER  IHPROTEHEHTS 


The  model  does  not  begin  tc  compute  a  track  until  the 
TMA  time  defined  by  the  user  at  simulation  initiation  is 
exceeded.  This  implies  that  all  sensors  and  operators  are 
equal,  which  is  not  realistic.  Furthermore,  this  does  not 
allow  for  accurate  information  to  be  used  at  time  of  detec¬ 
tion  unless  the  TMA  time  has  already  been  exceeded.  For 
instance,  a  passive  sonobuoy  dropped  in  front  of  a  contact, 
producing  a  CPA  (closest  point  cf  approach)  for  its  initial 
detection,  would  have  good  track  information  that  would  not 
be  utilized  by  the  mcdel;  for  only  information  after  the  TMA 
time  is  used  in  determining  the  FACTOR. 

This  TMA  initiation  time  was  included  in  the  model 
because  use  of  actual  target  information  provided  too  accu¬ 
rate  of  a  fix,  even  with  bad  sensor  information,  for  a 
player  to  experience  a  realistic  prosecution.  The  Kalman 
Filter  Mcdel  eliminates  the  need  for  this  waiting  time  since 
the  model  receives  information  as  an  operator  would  see  it 
(that  is,  appparent  bearing  resulting  from  apparent  posi¬ 
tion,  including  navigation  error,  and  sensor  bearing  error)  . 
Thus,  the  Kalman  Filter  Model  allows  use  of  all  sensor 
information  with  appropriate  errors  to  provide  realistic 
prosecution. 

The  TMA  model  attempts  to  simulate  a  changing  area  of 
probability  (AOP)  with  improved  solutions  taken  from  the 
table  as  SE  increases  and  true  tearing  changes.  The  problem 
of  using  true  information  instead  of  apparent  has  already 
been  discussed.  In  addition,  the  improved  AOP  is  heavily 
dependent  upon  the  drawing  of  a  random  number.  It  has  been 
observed  in  actual  trainers  that  the  AOP  fluctuates  as 
randomly  as  the  random  number  generator,  regardless  of 
sensor  information;  providing  confusing  information  to  the 
player.  Again,  the  Kalman  Filter  Model  eliminates  this 


0  <  SE 


2S 


.  IS 


.  05S  i 


m 

. 


R  =  Actual  Range 
S  =  Actual  Speed 


0 


1 

TMJ 

BEARING 

SIGNAL  CHANGE  IN 
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IN  dB  ( •  E) 
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II 
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GOOD 
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will  be  calculated  and  displayed.  Between  TMA  changes,  the 
displayed  range  is  updated  each  simulation  cycle  to  show  the 
range  of  the  target  estimated  from  the  TMA  course  and  speed. 

If  contact  is  lost  for  a  time  greater  than  the  user- 
defined  track  loss  time,  the  TMA  solution  will  no  longer  be 
displayed.  At  any  subsequent  redetection  of  the  same  target, 
a  new  solution  will  be  generated  after  the  appropriate  time 
interval  has  passed. 
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